IVP's: Longer Examples 



The fish population in a lake is not reproducing fast enough and the 
population is decaying exponentially with decay rate k. A program is 
started to stock the lake with fish. Three different scenarios are discussed 
below. 

Example 1. A program is started to stock the lake with fish at a constant 
rate of r units of fish /year. Unfortunately after 1/2 year the funding is cut 
and the program ends. Model this situation and solve the resulting DE for 
the fish population as a function of time. 

Solution. Let x(t) be the fish population and let A = x(0~) be the initial 
population. Exponential decay means the population is modeled by 

x + kx=f(t), x(0~) = A (1) 

where /(f) is the rate fish are being added to the lake. In this case 



f(t) = 



for < f < 1/2 
for 1/2 < f. 



First, write / in 'u-format': /(f) = r(l — u{t — 1/2)). 
Next, take the Laplace transform and solve for X(s). 

F(s)=£(f)(s) = r -- r -e- s/2 . 

sX-x(0-)+kX = F(s) (s+k)X- A = T -(\-e- s/1 ) 

=* X(s) = ^-+ . T ; , (l-g- s/2 ). 
w s + k s(s + k) y ' 

To find x(t) we temporarily ignore the factor of e~ s/2 and take Laplace in- 
verse of what's left, (using partial fractions). 



The f-translation formula says 

-s/2 



(%hs) = u i , - im i< 1 -* 



fc(^-l/2)^ 
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Putting it all together we get (in u and cases format). 

x(t) = Ae- kt + r -(l-e- kt )-u{t-l/2) r -(l-e- k ( t - 1/ V) 

( Ae- kt + {{1 - e- kt ) for0<f<l/2 
\Ae~ kt - {{e- kt +g-*(*-i/2)) for 1/2 < t. 



Example 2. (Periodic on/off) The program is refunded and the have 
enough money to stock at a constant rate of r for the first half of each year. 
Find x(t) in this case. 

Solution. All that's changed from example 1 is the input function f(t). We 
write it in cases-format and translate that to w-format so we can take the 
Laplace transform. 



/(0 = 



r for0<£<l/2 

forl/2<f<l 

r for0<f<3/2 

for 3/2 < t < 2 



= r(\-u{t- l -) + u(t-\)-u{t- 3 -) + ...) 

The computations from here are essentially the same as in the previous ex- 
ample. 

£(f) = 5(i 
x 



e -s/2 +e -s 



e~ 3s/2 + ...) 



s+k 1 s(s+k) 

x(t) = Ae- kt + \ 

' Ae- kt 4 
Ae~ kt - 



x(t) = < 



(l-e- s/2 + e- s - ...) 
(1 - e~ kt ) - u(t - 1/2) (1 - e-^- 1 / 2 )) + . . . 

r„-kt 



( e -kt _ e -fc(t-l/2)) 



for < f < i 
for | < f < 1 



r r 

k k 



(e- kt - e -Kt-i/2) + . . . + g-fc(»-»)) for n < f < n + \ 



: ( e -kt _ e -k(t-i/2) + _ _ _ _ e -fc(f-«-i/2) ) for n + i < f < n + 1 
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Factoring out e gives: 

f Ae- kt + \ 



x(t) = 



r„-kt 
k 



er»(l 



Ae 



-kt 



r„-kt 



(1 



Jc/2 



pk/2 



+ e K 



e k _ e 3k/2 + _ _ _ + e nk) torn <t <n + l/2 



for n + 1/2 < t < n + 1. 



Note that the constant term r/k is only present during periods of stocking. 

Example 3. (Impulse train) The answer to the previous example is a little 
hard to read. We know from experience that impulsive input usually leads 
to simpler output. In this scenario suppose that once a year r/2 units of 
fish are dumped all at once into the lake. Find x(t) in this case. 

Solution. Once again, all that's changed from example 1 is the input func- 
tion f(t). The IVP is still given by equation (1). 

fit) = T - (5(t) + 3(t - 1) + 5(t - 2 )+S(t-3) + .. .). 

This is called an impulse train. Its Laplace transform is easy to find. 



F s) 



(1 + e- s + e- ls + e- is + . . .) 



3s 



One nice thing about delta functions is that they don't introduce any new 
terms into the partial fractions part of the problem. 

r 



sX(s)-x(0-)+kX(s) = 

=> X(s) = 
Laplace inverse is easy: 

c- 1 



(1 



,-2s 



+ 



,-3s 



+ ...). 



(1 + e- s + e- 2s + e- is + . . .) 



s + k 2(s + k) 



3s 



1 

s + k 



Thus, 

x(t) = Ae- kt + r -( 



= e 



kt 



s + k 



= u(t — n)e 



+ Lu( t _ i) g -fc(t-D + r -u(t- 2)e- k ^ 



-k(t-n) 



-u(t — 3)e 



-k(t-3) 



Here are graphs of the solutions to examples 2 and 3 (with A = 0, k = 1, 
r = 2). Notice how they settle down to periodic behavior. 



t 





1 2 3 4 1 2 3 4 

Fig. 1. Graphs from example 2 (left) and example 3 (right). 
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